Cooling Nanomechanical Resonator by Periodic Coupling to Cooper Pair Box 
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We propose and study an active cooling mechanism for the nanomechanical resonator (NAMR) 
based on periodical coupling to a Cooper pair box (CPB), which is implemented by a designed series 
of magnetic flux pluses threading through the CPB. When the initial phonon number of the NAMR 
is not too large, this cooling protocol is efficient in decreasing the phonon number by two to three 
orders of magnitude. Our proposal is theoretically universal in cooling various boson systems of 
single mode. It can be specifically generalized to prepare the nonclassical state of the NAMR. 



PACS numbers: 85.85.+j, 85.25.Cp, 42.50.Pq, 32.80.Pj 
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Introduction — Recently nanomechanical resonators 
(NAMR) have been fabricated with high quality factors 
from (10 2 to 10 5 ) and large fundamental frequencies (in 
the range of MHz -GHz) JJ HQ. The NAMR has been 
shown as a good candidate for exploring various meso- 
scopic quantum phenomena at the boundary between 
classical and quantum realms. Until now NAMRs have 
been used in generating entangled states 0] , demonstrat- 
ing quantum non-demolition measurement Q and pro- 
gressive quantum decoherence Q, and implementing a 
two qubit quantum gate Q ■ 

The quantum nature of NAMR has been exhibited by 
the accurate measurement near the standard quantum 
limit HOI- But in most cases, to fully utilize the quan- 
tum wealths provided by NAMR it is necessary to cool 
the NAMR to its ground state. There have been some 
schemes for cooling NAMR [HI [H E3 Some of them 
are based on coupling with Josephson junction (J J) qubit 
and make use of feedback control and sideband cooling 
techniques. Characterized by the maximal ratio between 
the average number of phonons before and after cooling, 
the highest efficiency of some schemes 113 can be 
achieved when the initial number of phonons Nth of the 
NAMR is large enough (N th > 1CT 1 ). The cooling effect 
is evidently decreased as N t h gets very small. 

Motivated by the existing investigations mentioned 
above, we suggest a straightforward mechanism to cool 
the NAMR. Our scheme is also based on the coupling 
with the Cooper pair box (CPB), which is considered as 
a controllable two-level system. Different from the ex- 
isting schemes, our cooling protocol works efficiently in 
small Nth regime (Nth S 10 1 ). In our scenario, the in- 
teraction takes place periodically between the CPB and 
the NAMR. Before the interaction takes place in each 
cycle, the CPB is always set to its ground state so that it 
can absorb some energy from the NAMR during interac- 
tion period. A similar method has been used by to cool 
the microwave cavity [lii. Il4j. In principle, the present 
proposal can be generalized for cooling any single mode 
boson system via the coupling with a two-level system. 

One can intuitively compare our cooling mechanism 
with a classical analog. To cool a thermal box one can 



put a piece of ice into it and then drain the melt down 
water. It will take away part of the heat in the box. 
Naively, one can freeze the drained water into ice outside 
the box in some way and then place it back into the 
box. Repeat this process again and again until the box 
reaches the desired temperature. In our scheme, the CPB 
prepared in the ground (excited) state can be imagined as 
the ice (water) and the NAMR as the thermal box in the 
classical analog. However, due to quantum coherence, 
the mechanism of our cooling protocol is not as naive as 
this "ice-and-box" analog, because the energy loss of box 
is irreversible due to the second law of thermodynamics. 
The substantial difference of our protocol from the above 
classical analogue is the coherent oscillation of energy 
exchange between the qubit and the bosonic mode. In a 
sense, our scheme is more related to a type of quantum 
heat engine 




FIG. 1: The coupling system of the NAMR (within the rect- 
angle of dashed lines) and the CPB: the bang-bang coupling 
is implemented by a predesigned periodic series of magnetic 
flux <& x threading through the CPB. 

Model for our cooling protocol - - In our NAMR- 
CPB composite system shown schematically in Fig.l, 
the NAMR is directly connected to a CPB consisting of 
two Josephson junctions. The external magnetic flux <fr x 
threads through the SQUID, which can be used to ad- 
just its effective Josephson energy. C,j is the Josephson 
capacitance, V g the tunable gate voltage, C g the gate 
capacitance, V x the bias voltage on the resonator, and 
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C the effective capacitance of the NAMR. C x denotes 
the distribution capacitance between the CPB and the 
NAMR. Let E c = e 2 /(2 (C x + C g + Cj) be the charge 
energy unit and n g (x) = (C (x) V x + C g V g ) /2e be the 
total gate charge. 

The dependence of charging energy AE C (n c — n g (x)) 2 
on x results in the coupling between the NAMR 
with free Hamiltonian p 2 /(2m) + muj\x 2 /2 and the 
CPB with controllable Josephson tunnelling energy 
— 2Ej cos (7r$2;/$o) cos 9. Here, n c denotes the number 
of the excess Cooper pair on the island while its con- 
jugate variable is the phase difference 9 of the two 
sides of each junction. Usually the fluctuation of x is 
much smaller than the distance d between the NAMR 
and the CPB. At the exact resonance point defined by 
(C X V X + C g V g ) /2e = N/2 where N is odd, the CPB acts 
as a two level system, a qubit. We denote two linearly in- 
dependent charge states by |1) = \n c = (N — l)/2) and 
1 2) = \n c = (N — l)/2). Under the two level approxi- 
mation and the rotation wave approximation (RWA), we 
write the above Hamiltonian as the Jaynes-Cummings 
(JC) form H3 

H = E j cos ( ir-jr- J &z + u)Q(Ja + g (aa + + h.c) . (1) 

Here, cr z (er+,er_) are defined with respect to the new 
basis {|e) = (|1) + |2)) /y/2, \g) = (|1> - |2» /V2}; a (at) 
are the phononic creation and annihilation operators of 
the NAMR mode with the effective coupling constant 
g = 4E c n x yl/ 2mwo / d where n x = C x V x /(2e). 

In our protocol, the above JC type interaction is as- 
sumed to take place periodically with rate r a . Then it 
is switched off after the duty cycle interval r of the or- 
der of g^ 1 . This on-and-off switching can be realized by 
the magnetic flux (f> x . In fact, at the exact resonance 
point, the JJ tunnelling energy Ej cos (ir$ x /$o) is just 
the energy level spacing. It will be switched to the value 
resonant with the NAMR during the duty cycle r and 
far-off resonant outside this period. The similar manip- 
ulation has been used to create the non-classical photon 
state based on the superconductor devices 

To implement our proposal we have two main tasks to 
realize the periodical coupling: (a). Switching on and off 
the interaction between the CPB and the NAMR peri- 
odically, and (b) . Preparing the CPB to its ground state 
before the interaction takes place in each cycle. Both 
tasks can be accomplished via tuning the magnetic flux 
$ x with time. In fact, the gate charge fluctuation induced 
relaxation rate V (cj) = Tta n uj [coth (tu/2kfiT) + 1] /2 of 
the CPB at temperature T |l2j can be well controlled 
by varying the flux $ x since w = Ej cos (7r$ x /$o)- 
Here, ks is the Boltzmann constant, a g is about 
2e 2 R (C 2 + Cf) j (tt(C x + C g + C.j) 2 ) , and R the fluc- 
tuation impedance of V g and V x . Outside the duty cycle 
r, one can switch the energy spacing to of CPB to a large 



value to satisfy u> ± ujo >> g. In this case, the CPB- 
NAMR interaction is effectively switched off because of 
the far-off resonance (this can be deduced without RWA) 
while the decay process is enhanced. Thus the CPB can 
be prepared well in its ground state \g) for the up-coming 
interaction period. 

Master equation approach in steady states — We 
assume the coupling strength g is much stronger than 
r (u> a ) and k where K is the decay rate of NAMR, and 
the interaction period r is so short that T(cu a )T << 1, 
kt << 1. In this case both the decay of CPB during 
the duty cycle and the NAMR-environmcnt cou- 
pling can be omitted. Therefore, if the interaction is 
switched on at instance U, the reduced density operator 
p (ti + t) of the NAMR after a time interval r can 
be obtained through the action of the superoperator 
M (t) on the reduced density operator p(ti) at instance 
ti, i.e. p(ti+r) = M (r) [p (ti)], which is defined as 

M (r) [p (t,)] =Tr a [expH^V (tj) ® \g) {g\ exp(^r)] Tr a 

denotes tracing over the variables of CPB. h = gaa + + h.c 
is the JC type Hamiltonian at resonance in the in- 
teraction picture. Without any dissipation, the exact 
solution of the resonant JC model gives the explicit 
recursions 

Pn {ti +t) = |c g ,„ (t)| 2 p n (U) + |c ej „ (t)| 2 p n+ i (U) (2) 

for the diagonal elements p n = (n\ p \n) of p {U + r) for 
the number \n) state of NAMR phonon. Here, c e .„ (r) = 
sin (gT\/n + l) and c ffi „ (r) = cos (gT^/n) come from the 
exact solution of the resonant JC model. 

With the presence of the dissipation of NAMR, the 
evolution of p can be depicted by the course gained mas- 
ter equation 



^ = r a [M(r)-l}p + L[p]. 



(3) 



The super operator L in the above equation 
attributes to the dissipation and is defined 
as L [p] = —{n/2)Nt\ i (aa^p—2a'*pa-\-paa^ — 
(k/2) (JVth + 1) (a^ap — 2apa) + pa^aj where 
Ath = [exp (ujq /kgT) — 1] 1 is the average number 
of phonons of NAMR at temperature T before cooling. 
The master equation was initially presented for the 
case that the interaction between the two level system 
and the single mode oscillator is "turned on" randomly 
[HEi. In the case where the interaction in our scheme 
is periodically " turned on" , this equation can also lead 
to a correct stable solution. 

Without detailed computations, the average number of 
phonons (n) s =Tr[p s a^a] of NAMR in the steady state 
p s can be obtained from the above master equation: 

(n), = JVth - (r /«) An. (4) 

where An =Tr{a^a [1 — M (t)] p s }. Since 
Tr{a^a [1 — M (t)] p s } > can be proved with the 
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definition of M(r), we can conclude that (n) s < Nth, 
i.e., the NAMR can always be cooled when the steady 
state is reached. 

Dynamic process of cooling and the fluctuation of num- 
ber of phonons — The steady state solution of Eq. @ 
gives the phonon population 



N th l 



Pt, = Pn 17 



(5) 



Here, p!^ is n-th diagonal element of p s and Pq is de- 
termined by the normalization condition y^H n pf = 1 
|l7j . By a virtue of a numerical computations, we can 
describe exactly the evolution of vacuum state probabil- 
ity po and average number (n) of phonons in detail. Fig. 2 
demonstrates these results with experimentally rational 
parameters Nth ~ 1.7, r a /n = 133 and gr = n/8. It 
is shown that, at the time t a satisfying r a t a ~ 60, the 
steady solution JSJ is reached. In this sense the distribu- 
tion of the number of phonons in the data does not vary 
significantly after time t a . 



average number (n) s of phonons after cooling is drawn 
against A t h in Fig. 3. It shows that, with r a /« = 10 2 
and |c e ,o| = lj we have (n) s rts 10~ 2 iV tll for iV tll < 1; 
with r a jn = 10 3 and |c 6) o| = 1, we have (n) s rts 10 _3 iV tll 
for 7V th < 10 and (n) s w 10- 2 7V th for N th < 10 2 respec- 
tively. 
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FIG. 2: The time evolution of the average phonon number 
(n) and the vacuum state probability po. The unit of time is 
l/r Q . The fluctuation of phonon number in the steady state 
is shown in the inset. 

Apparently, we can improve the cooling effect by in- 
creasing r a . For a given value of r a / n, our cooling scheme 
works well when the average number of phonons N t h be- 
fore cooling is small enough, i.e., Nth << T a /n- Un- 
der this condition, we can achieve the maximum cool- 
ing effect by setting the duty cycle r = tt/ (2g) so that 
l c e,o| = 1- In this case the average number of phonons 
after cooling is 



(n) s «Pi« N t hK/r a . 



(6) 



This results implies that the number of phonons is re- 
duced by a factor of (r a /hi). However, if Nth is com- 
parable with (or larger than) |c e .o| 2 r a /n, our scheme 
does not work well. For JVth + 1 >> r a /K, we have 
Pn ~ Po [Nth/ (Nth + 1)]™, which implies that the aver- 
age number of phonons (n) s in the steady state is very 
close to Nth, the number of phonons before cooling. The 



FIG. 3: (color online) The cooling effect diagram. Here we 
assume a g m 1 x 10 -4 . The finial phonon number {n) s is 
drawn as a function of the initial phonon number N t h when 
Ta/n = 10 2 or 10 3 . In the solid (black) line, the thermal 



excitation of the CPB is not considered. In the dashed (blue) 
line, we assume the thermal excitation probability of the CPB 
is p — 10~ 4 . In the dotted (red) line, we assume p — 10 -5 . 

It is also noted that, since the interaction between the 
CPB and the NAMR is switched on and off again and 
again, and the master equation (j2J) is obtained via a 
coarse granulation approach, the number of phonons will 
have a fluctuation An =TrN {a^a [1 — M (r)] p s } even in 
the steady state. In the inset of Fig. 2, the fluctua- 
tion An is shown in the evolution curve of the number of 
phonons in the steady state. 

Experimentally feasible predications — We consider a 
NAMR with frequency lu =2tt x 10 2 MHz (0.5^eV) and 
quality factor Q = 2 x 10 5 . Our cooling protocol can 
be realized with the following experimentally accessible 
parameters, the decay rate k = irx 10 _3 MHz, the Joseph- 
son energy Ej w 4-7T x 10 4 MHz (lOOucV), the Coulomb 
charging energy Ec ~ 320/icV (i.e. ~ 250aF) ,the 
mutual capacitance C x = C g ~ 20aF, the impedances 
R w 50fi and the gate voltage V x w 0.25V. With these 
parameters, the interaction strength g is estimated to be 
2ir x 10MHz, the number of Cooper pair n x is about 15 
and Oi g m 1 x 10~ 4 . 

When the magnetic flux & x in the Hamiltonian JTJ 
is tuned to about 0.498<I>o on resonance, the interac- 
tion takes place. We set the interaction duration r = 
2.5 x 10~ 8 s so that |c e! o(r)| = 1. If the temperature 
T is 0.01 Kelvin, during the duty cycle, the decay rate 
r(w ) of CPB is about 0.56MHz. Just after the inter- 
action, the magnetic flux should be turned off to max- 
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imize the energy spacing. Then the CPB will decay to 
its ground state by a rather large decay rate T(Ej) m 40 
MHz. The CPB then decays quickly in r' = 0.275^ s. 
This elementary procedure is repeated at a frequency 
r a « 3 MHz. Namely, the interaction takes place ev- 
ery T\ = l/r Q = 0.3/xs. The time is sufficiently long for 
completing the interaction and the relaxation in each el- 
ementary procedure. Therefore, omitting the dissipation 
effect of the CPB, we can estimate from Eq. (JHJ) that the 
number of phonons after cooling has the magnitude of 
(n) s ~ N th /(r a / K )~W- 3 N th . 

In the above discussion, the parameter of crucial im- 
portance k is determined by the quality factor Q and 
the resonance frequency (Jo of the NAMR. As illustrated 
above, for Q ~ 2 x 10 4 , we end up with r a /n = 10 2 , 
and the cooling effect would be much worse. A care- 
ful analysis of the cooling cycle reveals that the cool- 
ing effect is also influenced by the dissipation of the 
CPB during the duty cycle. For high initial tempera- 
tures, we can replace the factor |c e! ;_i| 2 in formula (J2J) 
by |c e j-i| 2 F;_i where i^-i is the fidelity defined as 
= 1 — T (ujq) Jq \c e j-i (t')\ dt' from a first order 
perturbation calculation. This leads to a more precise 
phonon distribution p n in the steady state. For param- 
eters corresponding to a g « 1 x 10~ 4 as given above, 
the correction to p n due to the above perturbation ef- 
fect of the inclusion of CPB dissipation can be omit- 
ted when JVth S 10 2 - O n the other hand, the coupling 
of the CPB to the environment may cause the CPB to 
jump from the ground state to the excited state with a 
rate T (u>o), and subsequently emit a phonon via coherent 
Rabi oscillation. In the low temperature case, this proba- 
bility can be estimated as N t hT (u>o) J Q r |C e ,o (t — t')\ 2 dt' 
= A t hr (u>o) t/2. In our example as mentioned above, 
this value is about 7 x 10 _3 A t h- The detailed study on 
the influence of the dissipation effect of the CPB on our 
cooling protocol will be given in our future work. 

More generally, for the case of finite temperatures, 
there is always a thermal excitation for the CPB. In fact, 
even when off duty, the CPB decays to a mixed state 
with probability p = 1/ (1 + cxp[Ej/ksT]) in the ex- 
cited state |e) and (1 — p) in the ground state \g). For a 
finite value of p, the elements of the steady state density 
matrix p s should be modified to be 



Po 



N th l +p\c e ,i-i\ r a /K 
fJl {Nth + 1)1 + (1 -p) \ce,i-i\ 2 r a / K 



n 



(7) 



this leads to a reduction of the the cooling efficiency. In 
Fig. 3, the influence of this thermal excitation probability 
is illustrated. In the low temperature limit, the average 
number of phonons emitted by the CPB via this mech- 
anism during the duty cycle is still p. Combined with 
the analysis of the previous paragraphs, the lower limit 
phonon number fluctuation at steady state due to ther- 
mal excitations and dissipations of the CPB can be es- 



timated as 1/ (I + cxp[Ej/k B T}) + N th T (oj ) t/2, which 
vanishes in the limit of large Ej and low T(u>q). This 
can also be considered as the lower limit of the average 
phonon number after cooling. 

Conclusion with remarks about relations to maser - 
We should notice that the periodic cooling of the NAMR 
can be understood as an "inverse" of the "maser" mech- 
anism. In the usual maser process, the input excited 
atoms (molecules) can coherently heat the cavity and 
then coherently accumulates photons in a single state to 
enhance cavity field in a quantum way. In the present 
protocol, the NAMR is cooled by the CPB in its ground 
state. It should be emphasized that, with the same setup 
and operations similar to the above protocol, a "NAMR 
maser" can be devised if the CPB is initially prepared in 
its excited state |e) before the interaction is switched on. 
Namely, we can prepare the NAMR in the non-classical 
state with the number of phonons in super-Poissonian 
or sub-Poissonian distribution. On other hand, in many 
protocols of two-qubit quantum logic gates based on dif- 
ferent physical systems, the bosonic mode in its ground 
state can serve as a quantum data bus to transfer quan- 
tum information from its coupled qubit to another, or to 
entangle two qubit at a distance. Therefore, the protocol 
in this paper may have potential applications in quantum 
information theory. 
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